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Abstract 

We discuss a symmetry-adapted algebraic (or vibron) model for molecular spectroscopy. The 
model is formulated in terms of tensor operators under the molecular point group. In this way, we 
have identified interactions that are absent in previous versions of the vibron model, in which the 
Hamiltonian is expressed in terms of Casimir operators and their products. The inclusion of these 
new interactions leads to reliable spectroscopic predictions. As an example we study the vibrational 
excitations of the methane molecule, and compare our results with those obtained in other algebraic 
models. 
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1 Introduction 



The development and refinement of experimental techniques in high resolution spectroscopy has generated 
a wealth of new data on rovibrational spectra of polyatomic molecules. Highly symmetric molecules, such 
as tetrahedral XY4 systems, form an ideal testing ground. On the one hand, the high degree of symmetry 
tends to reduce the complexity of the spectrum and on the other hand, the use of symmetry concepts 
and group theoretical techniques may help to interpret the data and eventually suggest new experiments 
jj], I). A good example is provided by the methane molecule, for which there exists a large amount of 
information on vibrational energies. 

Ab initio calculations for rovibrational spectra of molecular systems attempt exact solutions of the 
Schrodinger equation. These calculations involve several configurations associated with the molecular 
electronic states and yield the force field constants Q from which the spectrum can be generated j^j. 
For small molecules this procedure is still feasible, but this is in general not the case for polyatomic 
molecules, due to the large size of the configuration space. Despite the progress made in ab initio 
calculations, a direct comparison with experimental vibrational energies of methane still shows large 
deviations, especially for vibrational states with a higher number of quanta. 

An alternative method is provided by algebraic (or vibron) models (for a review see || [?]]). The 
general method consists of two ingredients: (i) the introduction of U(k + 1) as the spectrum generating 
algebra for k degrees of freedom, and (ii) for a system of bosonic degrees of freedom the states are assigned 
to the symmetric representation [N] of U(k + 1). 

In its original formulation ]8|, 0], rotations and vibrations were treated simultaneously in terms of 
coupled U(A) algebras: Q = U\ (4) ® U2 (4) <S) ■ ■ ■ , by introducing a t/(4) algebra for each bond {k = 3). 
The electronic degrees of freedom can be included by introducing a unitary group for the electrons fl(i[| . 
For polyatomic molecules it was found to be more convenient to first separate the rotations and vibrations 
and subsequently to treat the vibrations in terms of coupled U(2) algebras [T^|: Q = C/i(2)<g)t/2(2)<g>. . . , 
introducing a U(2) algebra for each interatomic potential (k = 1). In this version of the vibron model 
the calculation of matrix elements is greatly simplified. An additional advantage is that it is well-suited 
to incorporate the underlying discrete symmetries. 

In a different approach, it has been suggested to use a U(k + 1) model for the k = 3n — 3 rotational 
and vibrational degrees of freedom of a n-atomic molecule |l3] ]. This model has the advantage that 
it incorporates all rotations and vibrations and takes into account the relevant point group symmetry. 
However, for larger molecules the number of possible interactions and the size of the Hamiltonian matrices 
increase rapidly. A similar approach can be used for the vibrations only pj . 

In this contribution, we discuss a symmetry-adapted version of the vibron model jy^ p^j [l8|| which 
is very well suited to describe the vibrations of polyatomic molecules, especially those with a high degree 
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of symmetry. The method is based on a set of coupled U(2) algebras, whose generators are projected on 
tensor operators under the molecular point group. In order to illustrate these ideas we first review the 
main ingredients of the U(2) vibron model, its connection with the Morse oscillator and the harmonic 
limit. Next we develop the formalism in more detail and take as an example the methane molecule which 
has tetrahcdral symmetry. Wherever possible, we make a comparison between the present formulation 
and other algebraic models. 



2 The U(2) vibron model 

The model is based on the isomorphism of the U(2) Lie algebra and the one-dimensional Morse oscillator, 



whose eigenstates can be associated with U(2) D 5*0(2) states |[iV],m) |19 . The U(2) D 5*0(2) algebra 
is generated by the set {G} = {N, J + , J_, J }, which satisfies the commutation relations 



Jo, J± 



±J± , 



2Jo , 



N,J, 







(1) 



with /i = ±,0. For the symmetric irreducible representation [N] of U{2), the Casimir operator is given 
by J 2 = N(N + 2)/4, from which follows the identification j = N/2. The SO (2) label is denoted by m. 



The Hamiltonian 



A 



corresponds to the Morse oscillator with energies 



= A ( f 2 J 2 



Em = A + 1) - m 2 ] 



AN 



(v 
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V 

iV 



(2) 



(3) 



where the label v = j — m denotes the number of quanta in the oscillator. The Morse eigenstates are 
denoted by \ [N], v) with v = 0, 1, .. ., [iV/2]. The first term in Em is the harmonic contribution, whereas 
the second term represents the anharmonicity which vanishes in the large N limit. 

The concept of the harmonic limit provides a link with a geometrical picture, and hence can be used 
to compare various models of molecular structure. Here we apply this procedure for the U(2) vibron 
model. The action of on the Morse eigenstates is 



J+\[N],v) = y/(N-v + l)v\[N},v-l) 



J-\[N],v) = y/(N-v){v + l)\[N],v + l) . 
2J \[N],v) = {N-2v)\[N],v) . 

Next we scale the generators J M appropriately and take the limit v/N — » 

J+ 



(4) 



lim ^\[N],v) 



ft\[N],v-l) = a\[N),v) 
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lim -^\[N},v) = VtT+T|[iV],« + l) = <f\[N],v) . 
}im n ^\[N],v) = (l-^)\[N],v) = (i-lja)\[N],v) , 

v/N— >0 iv iv iv 



iV 



|[JV],«) = |[JV],«) 



(5) 



In the harmonic limit (v/N — ► 0), the L/(2) algebra contracts to the Heisenberg-Weyl algebra which is 
generated by a^a, a\ a and 1. For example, we have 



1 

N 



*tl = 1 



(6) 



The other commutation relations can be treated similarly. The application of the harmonic limit of the 
Morse Hamiltonian of Eq. (0) gives 



lim -^-Hm 



lim 

iV— oo 2N 



J-J+ + J+J- 



A 1 

— (a^a + aa f ) = A(a) a + -) , 



with energies 



(7) 



(8) 



in agreement with the large N limit of Eq. (SI) 



3 Symmetry-adapted algebraic model 

The U(2) model described above was introduced to treat the stretching vibrations of diatomic molecules 
]TT| ]. For polyatomic molecules it was suggested to treat the vibrational excitations in terms of coupled 
U(2) algebras. This formulation was found to be very well suited to incorporate the underlying discrete 
symmetries ]l2[ . In particular, invariant interactions under the point group were constructed by applying 
proyection techniques on an expansion of the Hamiltonian in terms of Casimir invariants. In this section, 
we apply this process of symmetry adaptation to the generators of the U(2) algebras themselves, rather 
than to the Casimir operators. This procedure leads to new interaction terms. 

We illustrate the method by an application to the stretching and bending vibrations of methane. In 
the present approach, we associate a U(2) algebra with each relevant interatomic interaction. For the CH4 
molecule we have four U(2) algebras corresponding to the C-H interactions and six more representing 
the H-H couplings. The molecular dynamical group is then given by the product Q = U\(2) ® . . . 
t/io(2), where each Ui(2) algebra is generated by the set {Gi} = {-/V,;, J+,i, J-,i, Jo,i}, which satisfies the 
commutation relations 

[J Q ,i,J±,i] = ±J±,i, [J+,iJ-,i] = 2Jcm , [Ni,Jn,i] =0, (9) 
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with [i = ±,0. The labeling is such that i — 1, ... ,4 correspond to the C-H couplings while the other 
values of i are associated with H-H interactions ^0|. Here Ni is the i-th number operator. All physical 
operators are expressed in terms of the generators {G^}, and hence commute with the number operators 
iVj. For the CH4 molecule there are two different boson numbers, N s for the C-H couplings and Nb for 
the H-H couplings, which correspond to the stretching and bending modes, respectively. 

The tetrahedral symmetry of methane is taken into account by projecting the local operators {Gi}, 
which act on bond i, on the irreducible representations T of the tetrahedral group 7^. For the J^ i 
generators of Eq. (^|) we obtain the 7^ tensors 

10 

f l:i = ZXrA. (io) 

i=l 

where 7 denotes the component of T, and the label x refers to stretching (s) or bending (6). The 
expansion coefficients are the same as those given in [ p0[ for the one-phonon wave functions. The algebraic 
Hamiltonian is now constructed by repeated couplings of these tensors to a total symmetry A\. 

The methane molecule has nine vibrational degrees of freedom. Four of them correspond to the 
fundamental stretching modes {A\ © F2) and the other five to the fundamental bending modes (E © F2) 
pl| . A convenient labeling for the vibrational levels of CH4 is provided by (1/1,1/3 ,2/3,2/4 ), where v\, 
V2, vz and 1/4 denote the number of quanta in the A\ >a , Eb, F% tS and i<2,6 modes, respectively. The labels 
U are related to the vibrational angular momentum associated with degenerate vibrations. The allowed 
values are I4 = Ui, v\ — 2, . . . , 1 or for z/j odd or even The projected tensors of Eq. (|lo| ) correspond 
to ten degrees of freedom, four of which (Ay © F2) are related to stretching modes and six (Ay © E © F2) 
to the bendings. Consequently we can identify the tensor ^ as the operator associated to a spurious 
mode. This identification makes it possible to eliminate the spurious states exactly. This is achieved by (i) 
ignoring the T^\ b tensor in the construction of the Hamiltonian, and (ii) diagonalizing this Hamiltonian 
in a symmetry- adapted basis from which the spurious mode has been removed following the procedure 

of nn. 

3.1 Zeroth order Hamiltonian 

According to the above procedure, we now construct the Td invariant interactions that are at most 
quadratic in the generators and conserve the total number of quanta 

x 7 

Vr, = ^E^;r r ;. (11) 

X 7 
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Here T = A±, Fi for the stretching vibrations [x — s) and T = E, F% for the bending vibrations (x = b). 
In addition to Eq. (pi]), there are two stretching-bending interactions 

*- = ^IX^fi "- M 

The zeroth order vibrational Hamiltonian is now written as 

Ha = Wi Hai, s + w 2 + ^3 Kf 2 s + w 4 7Y_f 2 b + u 3 4 H sb 

+a 2 V Eb + a 3 V F2 s + a 4 Vf 2 ib + "34 V sb . (13) 

The interaction Va 1 s has not been included, since the linear combination 

i 4 

r=Ai,F 2 s i=l 

corresponds to the constant contribution N s + 2 to the energies. Similarly, for the bending vibrations the 
sum of the terms 

1 10 

E C^+VrJ = 7r r E A W + 2 )' ( 15 ) 

r=Ai,B,F 2 6 f=S 

corresponds to a constant 3(iVb+2)/2. However, in this case the interactions Hai b an d Va 1 b have already 
been excluded in order to remove the spurious contributions from the Hamiltonian. 

The Hamiltion of Eq. (|l^) is equivalent to an expansion in terms of Casimir operators. It has the 
advantage, though, that the spurious contributions have been eliminated from the outset. A comparison 
with the Hamiltonian of (2Cj yields three conditions on their parameters 

A 5 + 2_B 5 ,io + 8B 5 6 = , 

Bl,5 + -Bl,8 = , 

Ai,6 + Ai l8 = 0. (16) 

The first condition eliminates the spurious interaction from the bending Hamiltonian of p0| , whereas the 
latter two eliminate the spurious contributions from the stretching-bending interactions. We note that the 
condition on the Hamiltonian that was used in p0| to exclude the spurious terms, does not automatically 
hold for states with higher number of quanta, nor does it remove all spurious contributions. 

3.2 Harmonic limit 

In the harmonic limit the interaction terms of Eq. ([l3]) have a particularly simple form, which can be 
directly related to configuration space interactions 

^tfr. = 5EK ,tfl ? ,+fl ?*< ,t ) ' 

7 
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lim Vr T = , 

N x — >oo 



lim H sb = I^( a ^t a ^, +a F 2iSa F 2 , 6t 



N a .N b ->co 2 

7 

lim V sb = . (17) 

N s ,N b -*oo 

Here the operators aE* are given in terms of the local boson operators a% through the coefficients a " 
given in Eq. ( |l0|) 
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«7 X = E a 7>*> ( 18 ) 
i=l 

with a similar relation for the creation operators. From Eq. (|l7]) the physical interpretation of the inter- 
actions is immediate. The 7ir x terms represent the anharmonic counterpart of the harmonic interactions, 
while the Vr^. terms are purely anharmonic whose contribution to the excitation energies vanishes in the 
harmonic limit. 

We note, that the recently introduced boson-realization model [^3| corresponds to the harmonic limit 
of the present approach, since it is formulated directly in terms of the boson creation and anniliation 
operators, aE"' and a^ x . The difference between the two lies in the anharmonic contributions which are 
implicit in the U(2) approach, but which vanish in the harmonic limit p2| . 

3.3 Higher order interactions 

The zeroth order Hamiltonian of Eq. (jl^) is not sufficient to obtain a high-quality fit of the vibrations 
of methane. For example, the results presented in pC| ] were obtained by fitting 19 vibrational energies 
with a r.m.s. deviation of 12.16 cm" 1 . The boson-realization model of which, as was shown above, 
corresponds to the harmonic limit of the present approach was applied to the same 19 vibrations with a 
r.m.s. deviation of 11.61 cm -1 . We note, however, that the latter calculation includes some higher order 
interactions, without significantly improving the results. 

Several physically meaningful interaction terms that are essential for an improved fit are not present 
in Eq. (|l3|). They arise in the present model as higher order interactions. Products of Hi and Vj 

HiHj , ViVj , , (19) 

are equivalent to an expansion in powers of Casimir operators. These terms only involve intermediate 
couplings with T = Ai symmetry, since Hi and Vj themselves are scalars under the tetrahedral group. 
However, there exist other interaction terms that involve intermediate couplings with r = A2, Fi, E, F2 
symmetry. For example, the interactions 

522 l M l M + 533 £ f& f& + 544 J2 ^ ^7 + 534 E ^7 ft > ( 20 ) 
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with 

l M = -iV2— [f Eb x ff b ] M , 
Nb 

1* = +iV2^-[f^xf^ . (21) 

split levels with the same (yi, V3, but with different I2, h and/or I4. The square brackets in 
Eq. ( pl| ) denote the tensor couplings under the point group 7^. Similarly, all higher order terms and 
anharmonicitics can be constructed in a systematic way. Each one of the interaction terms has a direct 
physical interpretation and a specific action on the various modes. 

For the study of the vibrational excitations of methane we propose the following 7^ invariant Hamil- 
tonian jl§ |§ 

H = UJi HAi >s + HE b + ^3 Wjr 2 s + CO4, Hf 2 , b + a 3 Vf 3 , s 



-X23 ( 7~(-E b 7~Lf 2 . s ) + -^24 ( 7~(-E b Wf 2 s ) + X34 I Hf 2 _, Hir 2 



+522 l M l M + 533 E 'ft + ^4 E ft f £ + ^ E ^ ^7 

7 7 7 

+^33 Css + ^44 Cfob + ^34 Csb + ^23 &2S + ^24 026 ■ (22) 

The interpretation of the u>i and 0:3 terms follows from Eq. (|l7|). The Xy terms are quadratic in the 
operators 'Ht x and hence represent anharmonic vibrational interactions. The gij terms are related to 
the vibrational angular momenta associated with the degenerate vibrations. As mentioned above, these 
interactions, which are fundamental to describe molecular systems with a high degree of symmetry, are 
absent in previous versions of the vibron model in which the interaction terms are expressed in terms 
of Casimir operators and products thereof. In the harmonic limit, the expectation value of the diagonal 



terms in Eq. (22) leads to the familiar Dunham expansion J2l| 



E uh (vi + §) + E E *i («< + y)fe + y) + EE 9ij kh ■ (23) 

Here di is the degeneracy of the vibration. The terms in Eq. give rise to further splittings 
of the vibrational levels v 2l v 3 , V4) into its possible sublevels. In the harmonic limit the terms 
have the same interpretation as in B. The O ss , Obb and O s b terms give rise to a splitting of the E 
and F 2 vibrations belonging to the (ui, vlf, v l s 3 , v l £) = (0, 0°, 2 2 , 0°), (0, 0°, 0°, 2 2 ) and (0, 0°, l 1 , l 1 ) levels, 
respectively. Similarly, the C>2s and C>2b terms split the F\ and F2 vibrations belonging to the (0, l 1 , l 1 , 0°) 
and (0, l 1 , 0°, l 1 ) overtones, respectively. 



4 Results 



The Hamiltonian of Eq. (|22) involves 23 interaction strengths and the two boson numbers, N s and 
Nb- The vibron number associated with the stretching vibrations is determined from the spectroscopic 
constants uj e and x e to e for the CH molecule to be N s = 43 p5| . The vibron number for the bending 
vibrations, which are far more harmonic than the stretching vibrations, is taken to be = 150. We 
have carried out a least-square fit to the vibrational spectrum of methane including 44 energies. We find 
an overall fit to the observed levels with a r.m.s. deviation which is an order of magnitude better than in 
previous studies. While the r.m.s. deviations of (2(| and |2j| are 12.16 and 11.61 cm -1 for 19 energies, 
we find a r.m.s. of 1.16 cm" 1 for 44 energies. The values of the fitted parameters as well as all predicted 
levels up to V = 3 can be found in |ll| |24| . 

The «3 term plays an important role in the calculation. It is completely anharmonic in origin and its 
contribution to the excitation energies vanishes in the harmonic limit. In order to address the importance 
of this term in Eq. ( |22| ) we have carried out another calculation without this term. With one less 
interaction term the r.m.s. deviation increases from 1.16 to 4.48 cm" 1 . This shows the importance of 
the term proportional to 013 to obtain an accurate description of the anharmonicities that are present in 
the data. The absence of the 013 term in the second calculation can only partially be compensated by the 
anharmonicity constants Xij. 



5 Summary and conclusions 

In summary, we have discussed a symmetry-adapted algebraic model for molecular vibrations, in which 
the symmetry adaptation is applied at the level of the generators. This procedure has several interesting 
aspects: 

• it provides a systematic procedure to construct all interaction terms up to a certain order, 

• the harmonic limit gives a relation with configuration space interactions and Dunham expansions, 

• the spurious states can be removed exactly. 

The application to the 44 observed vibrational excitations of methane gives a good overall fit with a 
r.m.s. deviation of 1.16 cm -1 corresponding to an accuracy of ~ 0.01 — 0.10 %, which can be considered 
of spectroscopic quality. 

It was pointed out that the Vf 2 s term in combination with the anharmonic effects in the other 
interaction terms plays a crucial role in obtaining a fit of this quality. Purely anharmonic terms of 
this sort arise naturally in the symmetry-adapted algebraic model, but vanish in the harmonic limit. 
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Physically, these contributions arise from the anharmonic character of the interatomic interactions, and 
seem to play an important role when dealing with molecular anharmonicities. 

We have established an explicit relation with the algebraic model of , in which the Hamiltonian is 
expressed in terms of Casimir operators. A comparison between the two methods yields three constraints 
on the parameters, which remove the spurious components from the Hamiltonian of pOfl . A comparison 
with the boson-realization model of |2^| shows that this model corresponds to the harmonic limit of the 
present approach. 

The predictability has been tested by systematically adding levels with higher number of quanta in 
the fitting procedure. The slow variation in the parameters shows that the model has a high degree 
of predictability. The application to methane |l8| and to other molecules [H| [l7| suggest that the 
present model provides a numerically efficient tool to study molecular vibrations with high precision 
(r.m.s. deviations of ~ 1 cm -1 ). 
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